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We study the integrated Sachs-Wolfe (ISW) effect in ghost-free, massive bigravity. We focus on 
the infinite-branch bigravity (IBB) model which exhibits viable cosmic expansion histories and stable 
linear perturbations, while the cosmological constant is set to zero and the late-time accelerated 
expansion of the Universe is due solely to the gravitational interaction terms. The ISW contribution 
to the CMB auto-correlation power spectrum is predicted, as well as the cross-correlation between 
the CMB temperature anisotropies and large-scale structure. We use ISW amplitudes as inferred 
from the WMAP 9-year temperature data together with galaxy and AGN data provided by the WISE 
mission in order to compare the theoretical predictions to the observations. The ISW amplitudes in 
IBB are found to be larger than the corresponding ones in the standard ACDM model by roughly a 
factor of 1.5, but are still consistent with the observations. 
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I. INTRODUCTION 

The theoretical and phenomenological implications of 
massive gravity and its bimetric extension have been 
extensively studied since the ghost-free, nonlinear formu¬ 
lation was constructed in 2011 m , some 70 years after 
the linear theory of massive gravity was proposed by Fierz 
and Pauli m (see Ref. m for a comprehensive review). 
It has been shown that in order for the nonlinear theory to 
avoid the so-called Boulware-Deser (BD) ghost [12], one 
needs to extend the gravitational sector by introducing 
at least one extra metric-like spin-2 tensor field. This 
so-called reference metric can either be fixed, resulting in 
a theory of massive gravity where gravitons possess five 
degrees of freedom, or be given dynamics, giving rise to 
a bimetric theory with two interacting tensor fields and 
seven degrees of freedom, corresponding to one massless 
and one massive graviton. The two theories are often re¬ 
ferred to, respectively, as the de Rham-Gabadadze-Tolley 
theory of massive gravity, or dRGT, and the Hassan-Rosen 
theory of massive bigravity. 

Soon after it was proposed, the original dRGT theory, 
where only the dynamical metric couples to matter, was 
shown to suffer from a no-go theorem forbidding flat and 
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closed Friedmann-Lemaitre-Robertson-Walker (FLRW) 
cosmological solutions on a flat reference metric [T3]. So¬ 
lutions with an open FLRW metric or a non-flat reference 
metric do not save the theory, as they suffer from the 
so-called Higuchi instability [14] or other types of insta¬ 
bilities [T5H50] . Several solutions to these problems have 
been proposed in the literature by either adding extra de¬ 
grees of freedom to dRGT [QJ fSTUSTI (for example stable 
and self-accelerating cosmological solutions are possible in 
the quasidilaton extension of massive gravity [551 EH| ) or 
assuming inhomogenous and/or anisotropic metrics which 
are arbitrarily close to FLRW solutions on observable 
scales [HUI31E03S2]. As a different route to overcome 
the no-go theorem, it has recently been shown in Refs. E3- 
[55] that exact FLRW solutions in dRGT can exist if at 
least some matter couples to both the dynamical and 
the fixed reference metrics through a generalized matter 
coupling (see, however, Ref. [35] for complications this 
approach would introduce for cosmological studies). The 
hybrid metric introduced in this extended theory does, 
however, reintroduce the BD ghost above a certain energy 
scale []] which is a generic feature of theories where both 
metrics couple to matter simultaneously GEBEllSEHni, 
including the simple case of coupling the two metrics to 
matter minimally [321 - H5] . 


1 This scale is, however, thought to be above the strong coupling 

scale of the theory. Therefore phenomenology of the theory, 

including its cosmological implications, can likely be safely studied 

below a cut-off scale larger than the strong coupling scale ESI- 
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In this paper we consider the Hassan-Rosen theory of 
massive bigravity in its original formulation, i.e., where 
both metrics are dynamical while only one of them couples 
to matter. Many of the problems mentioned above are 
cured in bigravity, making it an attractive modification to 
general relativity (GR) with interesting phenomenology. 
The cosmology of bigravity has been extensively studied 
at both the background and perturbative levels (see, e.g., 
Refs. USUI]), and the theory has been shown to be con¬ 
sistent with observed lensing and dynamical properties of 
local sources such as elliptical galaxies [52]. In particular, 
it has been shown that the theory can provide viable cos¬ 
mological models as alternatives to the standard ACDM 
model, consistent with all existing cosmological observa¬ 
tions while offering clear signatures that can be verified 
using future large-scale structure experiments (see also 
Refs. [621464] for cosmological studies of bigravity in the 
context of the generalized matter coupling). Perturbative 
studies have placed strong conditions on the parameter 
space of the theory as a large fraction of the models 
generate early-time instabilities which, although they do 
not necessarily rule these models out, complicate their 
comparison to observations as linear perturbation theory 
cannot be employed. A specific submodel, infinite-branch 
bigravity (IBB) [59], exhibits both viable background so¬ 
lutions and stable linear perturbations that are consistent 
with observational data. This is the model that we focus 
on in this paper]^] 

As a continuation of the cosmological studies of bi¬ 
gravity, here we perform an analysis of the integrated 
Sachs-Wolfe (ISW) effect [5B] (see Ref. [67] for a recent 
review) in the framework of the original Hassan-Rosen 
theory. The ISW effect is the modification of the cos¬ 
mic microwave background (CMB) spectrum due to the 
time evolution of gravitational potentials generated by 
the large-scale structure of the Universe. When CMB 
photons, traveling between the surface of last scattering 
and our detectors, enter and then leave a gravitational 
potential well, they blueshift and redshift, respectively. 
If the gravitational potential does not evolve with time, 
which is the case in a flat, matter-dominated universe 
governed by GR, photons do not lose or gain energy when 
leaving gravitational wells, and the CMB spectrum is un¬ 
modified. If the Universe is not flat or matter-dominated, 
the gravitational potential does evolve with time, result¬ 
ing in a net change in the photon energy. Therefore, in 
GR, the ISW effect is closely linked to dark energy (or 
curvature). As dark energy starts to dominate and the 
cosmic expansion starts to accelerate, the potential wells 
become shallower with time and CMB photons gain a 


2 During the completion of this paper, an update to Ref. m 
appeared which studied the Higuchi bound in IBB and showed 
that it is not satisfied (see also Ref. ESI which appeared later). 
This can potentially be dangerous since instabilities might appear 
at higher-order perturbations. Whether that is the case or not 
needs to be studied by analyzing nonlinearities; we leave the 
investigation of this for future work. 


net blueshift. Since the acceleration of the expansion 
becomes important at low redsliifts, the ISW effect is a 
late-tinre phenomenon]^] Measurements of the ISW effect 
therefore provide a complementary method to investigate 
the properties of dark energy. A similar effect is expected 
if the acceleration is due not to dark energy but rather 
to large-scale modifications to GR, such as in bigravity, 
where dark energy is replaced by non-standard gravita¬ 
tional effects. However, since the gravitational potential 
in these theories behaves differently at late times than 
in GR, the ISW effect can be different, and measuring 
its properties can provide another observational tool to 
distinguish modified gravity from GR and dark energy 
models. 


The ISW effect is normally not measurable directly from 
the CMB spectrum [61], since it is subdominant compared 
to the ordinary Sachs-Wolfe effect caused by gravitational 
redshift at the time of decoupling. Therefore in order to 
measure the ISW effect, one has to cross-correlate the 
CMB temperature fluctuations with other tracers of the 
gravitational potential, such as the large-scale distribu¬ 
tion of galaxies and quasars [65] ■ This has motivated 
several analyses trying to detect the ISW signal using ei¬ 
ther the angular cross-power spectrum between the CMB 
temperature anisotropies and large-scale structure, or a 
similar correlation in pixel space using stacking techniques 
(see, e.g., Refs. [551 [704457] and references therein). The 
significance of detection has been around 3<r on average, 
with the strongest detection (~4.5cr) reported in Ref. [83j . 
Interestingly, many of these previous studies reported a 
cross-correlation signal which is systematically larger, by 
about ler-2cr, than the value predicted in the standard 
ACDM model (see, e.g., Refs. [55H5H1551155] 1. with some 
studies going even further and claiming an observed signal 
> 3(7 above the predicted standard model value [5511100] 
using stacking techniques. Note, however, that later stud¬ 
ies detected a cross-correlation signal consistent with the 
prediction of ACDM [55]. 


This paper is organized as follows. In Sec. [H]we present 
the Hassan-Rosen theory of bigravity, and in Sec. |III| 
we review its predictions for cosmic expansion histories, 
focusing on the specific case of IBB. The perturbative 
framework is briefly presented in Sec. IV and in Sec. [V] 
we discuss the ISW effect and how it can be approached 
in massive bigravity. We derive the theoretical ISW con¬ 
tributions to the CMB power spectrum, as well as to 
the cross-correlation between the CMB anisotropies and 
large-scale structure. Our numerical results are presented 
in Sec. m where they are compared to the observational 
measurements of the ISW effect. We finally discuss our 
findings and conclude in Sec. |VII| 


3 There also exists a subdominant early-time ISW effect due to 
the slight impact of the radiation component on the background 
expansion at the time of decoupling. We disregard this effect in 
this paper. 
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II. THE THEORY OF MASSIVE BIGRAVITY 

The Hassan-Rosen theory of ghost-free, massive bigrav¬ 
ity has the Lagrangian [8] 

Ml , - M2 ,- 

C = -det gR g - det f Rf 

4 

+ TO 4 y/~ det g /3 n e n (Vs” 1 /) + A/^ det 9^ m. 

n—0 

(1) 

where and I?/ are the Ricci scalars for the metrics 
and respectively. The functions e n are the ele¬ 
mentary symmetric polynomials of the eigenvalues of the 
matrix \/g~ 1 f, where the matrix square root is defined 
such that \/g~ 1 f\/g~ 1 f = g^ v f^v The forms of these 
polynomials are given in, e.g., Ref. [8j. The quantities 
/3„ (n = 0,1,2, 3,4) are free parameters of the theory to 
be constrained by observations. As pointed out earlier, 
in this paper we consider the so-called singly-coupled bi¬ 
gravity theory where the matter Lagrangian C m couples 
only to g lxv . 

The generalized Einstein equations for the two metrics 
are 

3 1 

G% + m 2 J2 (- 1 )" Pn9„>Xtu)u 

n =0 9 

(2) 
3 

GU + m 2 * (-1)” fa-nUxYfc )v (/FV) = 0, (3) 

n—0 

where the Einstein tensors corresponding to g gv and 
are defined as usual, = R 9 — \R g g^ with R 9 
the ^-metric Ricci tensor and similarly for Gj iv . and the 

matrices Y( n \ are functions of y/ s” 1 / which are defined, 
e.g., in Ref. 0. is the stress-energy tensor of matter 
defined in the usual way with respect to the physical 
metric, g gi/ . Here we have performed a rescaling of the 
parameters and fields so that M'j becomes a redundant 
parameter]^] This is possible when we couple matter only 
to g^. Finally, by combining the Bianchi identity with 
the field equations and the conservation of the stress- 
energy tensor, we arrive at the Bianchi constraint on the 
interaction (mass) terms, 

(W7) = o, (4) 

n= 0 

where \7 g is the covariant derivative operator with respect 
to 9iiv 


4 See, however, Ref. cnn, which appeared during the completion 
of this work, for caveats associated with this rescaling. 


III. COSMOLOGICAL BACKGROUND 

Following the literature on the cosmology of massive 
bigravity, and assuming the cosmological principles of 
homogeneity and isotropy, we assume here that the Uni¬ 
verse at the background level can be described by FLRW 
metrics for both g^ v and /^[j For reasons of brevity, 
we do not present the details of how to derive the back¬ 
ground solutions, and instead refer the reader to, e.g., 
Refs Braun EH EDM®. 

The background has three free metric functions: the 
g- metric scale factor a g: the /-metric lapse Nf, and the 
/-metric scale factor a/. The g-metric lapse N g can be 
chosen arbitrarily through a choice of time coordinate; 
we work with conformal time in this paper, i.e., we set 
N g = a g . The matter sector is described by the energy 
density p and the pressure p, which satisfy the usual 
continuity equation with respect to the metric g^, 

P + 3— (p + p) = 0. (5) 

a g 

The Bianchi constraint 0 yields two possible solution 
branches. The first branch places an algebraic constraint 
on the quantity r = af/a gi which becomes a constant 
with a value that depends on the parameters /?i, /3 2 and 
/3 3 . The second branch places a dynamical constraint 
on the /-metric lapse: Nf = a g a,f/a g . Focusing on the 
second (dynamical) branch, which gives interesting non- 
GR solutions, and defining 


B o (r) 

= 0o + 3/3i r + 3f3 2 r 2 + p 3 r 3 , 

(6) 

Bi(r) 

■aS 

III 

3 + 3(3 2 r 2 + 3/3 3 r 1 + /3 4 , 

(7) 

the field equations become 



3 H 2 g 


(8) 


P 

Ml 

II 

to 

to 

to 

o 

(9) 


where H g = a g /a 2 is the Hubble factor for the physical 
metric g gl ,. The first equation is the standard Friedmann 
equation modified by a mass term that depends on r, 
while the second equation determines r algebraically. 

As mentioned in Sec. |TJ in this paper we focus on the 
IBB model where only /3\ and /3 4 are non-vanishing. This 
model contains no explicit cosmological constant term, 
has stable linear perturbations, and is a viable alternative 
to ACDM in terms of background expansion and growth 
of structure [55]. In IBB, r starts off from a large value at 
early times and decreases with time to a finite value in the 
far future. The best-fit parameter values for this model 
from the comparison of its theoretical predictions to the 


5 See Ref. 11021 and references therein for other possible metrics in 

bimetric cosmology. 
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cosmic expansion history, as well as the linear growth of 
structure, are /?i = 0.48 and /?4 = 0.94 [58lf59 ; j^] where the 
parameters are measured in units of H q with H 0 = H g 0 
the Hubble rate today, and the mass parameter m 2 is 
absorbed into the parameters /3 n . These are the parameter 
values which we use in this paper when comparing the 
ISW predictions to observations. 


IV. COSMOLOGICAL PERTURBATIONS 

Cosmological perturbations in massive bigravity have 
been derived and studied in several previous works. They 
were first presented in Ref. 03, and the scalar pertur¬ 
bations were then further studied in Ref. J4B] in the su¬ 
perhorizon limit and in connection with primordial per¬ 
turbations. The perturbative equations were rederived 
in Ref. [H] and studied in the future de Sitter regime. 
In Ref. [53], the linear growth of cosmological structure 
in the quasi-static limit and for subhorizon modes was 
studied for a minimal bigravity model, where only the 
parameter j3\ is non-vanishing. It was shown that this 
minimal model, despite possessing a viable expansion his¬ 
tory, is unstable at early times on small scales, confirming 
the instabilities found previously in Ref. m- An exten¬ 
sive analysis of the theory for scalar perturbations and 
for various combinations of the parameters was then per¬ 
formed in Ref. [58], the quasi-static limit for subhorizon 
scales was studied, several modified gravity parameters 
were calculated, and deviations from GR predictions mea¬ 
surable by future large-scale structure experiments were 
presented. Ref. [59] then studied the problem of instabil¬ 
ities in further detail, and identified specific submodels 
and parameter combinations which avoid the instabilities; 
this introduced the IBB model that we consider in this 
paper. The problem of gradient instabilities has been 
further studied in Ref. m- Finally, Ref. [BO] presented a 
comprehensive analysis of the scalar, vector, and tensor 
cosmological perturbations, confirming the results found 
in the previous studies. Cosmological perturbations have 
also been studied in some extensions of the theory; see, 
e.g., Ref. [5B] for the case where the / metric is coupled 
to a second matter sector. 

In order to obtain the perturbation equations used in 
our calculation of the ISW effect for bigravity, as will be 
discussed in the next section, we follow Refs. [55HB0] and 
use the following ansatz for the perturbative metrics: 

ds 2 g = — a 2 (1 + Eg) dr] 2 + 2 a 2 diF g di]dx 1 

+ a 2 [(1 -I- A g ) 5ij + didjBg] dx l dx j , (10) 


Note, however, that background observations can only pick out a 
curve in the (3i~(34 space, and while perturbations can in principle 
break this degeneracy, present data do so poorly, cf. Fig. 3 in 
Ref. |59| . Therefore while the parameters quoted here fit the best, 
a range of other parameters along this best-fit curve fits nearly 
as well. 


ds 2 = — N 2 (1 + Ef) drf + 2NfajdiFfdr]dx l 

+ a 2 [(1 -I- Af ) Sij + didjBf] dx^dxE ( 11 ) 


The equations of motion for the perturbation variables 
and their numerical solutions have been given in the 
same references (see also Refs. gTHM [52] [37]). Since 
the equations are rather lengthy, we refer the reader to 
those papers for details and here we only present some 
key definitions and assumptions used in our ISW analysis. 

Forming the gauge-invariant variables ESP] 



( 12 ) 

(13) 

(14) 


where p m is the background energy density of matter and 
5p m is its perturbation, the equations of motion can be 
written as two second-order differential equations for the 
gravitational potentials and fk. Having solved these 
numerically (in the gauge F g = Ff = 0), we can obtain 
the time evolution of the gauge-invariant density contrast 
of the matter perturbations. 

<f>, T and S are the quantities needed when computing 
the ISW effect. We again emphasize that the model we 
study in this paper is IBB, therefore the variables and 
the corresponding equations of motion we compute are 
the ones specific to the IBB model. 

Before moving on to the discussion of the ISW effect 
in bigravity compared to that in ACDM, it is illustrative 
to see how the density contrast of matter perturbations 
differs between the two models on different scales. We 
therefore plot in Fig. [T| the ratio of the density contrasts 
today computed in ACDM (with = 0.3) and IBB 
(with /?! = 0.48 and (3^ = 0.94). We have used the same 
initial conditions in both models for the gravitational 
potentials at the time of decoupling. The figure shows 
that it is only for modes close to or larger than the horizon 
scale today that the predictions of the two models differ 
significantly, with IBB predicting much less superhorizon 
structure than ACDM. As we will see in the following 
sections, the amplitude of the ISW effect is large only 
at low multipoles (£ < 100); the multipoles 2 < £ < 10 
receive contributions from the modes 1 < k < 100, where 
the scales k are measured in units of Hq. 


V. THE ISW EFFECT 


The ISW effect is due to the time evolution of gravita¬ 
tional potentials $ and T between the epoch of recombi¬ 
nation (last scattering surface) and today, generating the 
following secondary anisotropies on the CMB sky: 


Vo 


= I (# + $) cfy. 

1 J ] 


ISW 


(15) 





one has 


5 



*-5-f f I '<*■*>■ 


(19) 


where XL m is the matter density parameter at present time. 
Plugging these relations and equations in Eqs. © and 
(18), one recovers the known expressions for the weight 


functions in GR. 

In massive bigravity, $ is no longer equal to 4> at all 
times, and furthermore, the time evolutions of the density 
contrast and potentials depend on k. It is possible to 
normalize the weight functions with respect to either the 
density contrast today, or the potentials at the time of 
decoupling. We explore both options below. For the 
power spectrum, we use the standard form 


FIG. 1. Ratio of the gauge-invariant density contrast in IBB 
today to the one in ACDM, as a function of scale, computed 
using the same initial conditions for the gravitational potentials 
at the time of decoupling. Scales k are in units where Ho = 1. 


As mentioned in Sec. [TJ, the ISW effect is measur¬ 
able through the cross-correlation between the CMB 
anisotropies and the large-scale structure. Therefore, in 
order to compare the predictions of bigravity to obser¬ 
vations, we need to compute this cross-correlation. The 
standard quantity to work with is the cross-correlation 
angular power spectrum, Cj s , since most of the ISW 
analyses are done in harmonic space. Following Ref. [55] , 
Cj 9 is given by 

cj 9 = cf 9 = 4tt f l(k)K}(k)K 9 (k). (16) 

This form for Cj 9 assumes that there are no other contri¬ 
butions to the correlation between the CMB temperature 
anisotropies and large-scale structure at low redshifts. 
Here, Kf and K 9 are weight functions given by 


A 2 — <f 2 

— °H 


Ho 


n+3 


r 2 (fc), 


( 20 ) 


where 5h is the amplitude of the density contrast at the 
horizon scale today, n is the spectral index, and T(k) is 
the transfer function. For the transfer function we use 
the Bardeen-Bond-Kaiser-Szalay (BBKS) approximation 
m- Normalizing the density contrast with respect to 
the potentials at the time of decoupling, one can use Eq. 
(19) to infer the spectrum of 4> at that time. One then 


replaces A^ with A| in Eq. (16), and normalizes the 


weight functions with respect to the potentials at the time 
of decoupling instead of the density contrast today. In 
GR, both approaches obviously give the same results, but 
in massive bigravity Cj 9 varies slightly depending on the 
choice of the normalization scheme. 

Besides the cross-correlation spectrum Cj 9 , one can 
also compute the ISW auto-correlation spectrum C] sw 
directly, although, as mentioned earlier, this quantity 
cannot be used to detect the ISW effect from observations. 
Cf w is given by 


C) sw = Cr = 4t r 




dk 

t A| (k)Kf(k)I<f(k). (21) 


K 9 (k) = 
Kf{k) = 


dzh i Z )~J~ [ k X( z )\ 


'dz± 

dz 


dz 6(k , 0 ) 

<I> + 4/ 


S(k,0) 


3t [kx(z)\, 


(17) 

(18) 


where b(z) is the bias factor which relates tracers to 
the underlying mass distribution, dN/dz is the redshift 
distribution of the density tracers normalized so that 
f dz^p- = 1, and jt. are the spherical Bessel functions 
with arguments depending on x{ z ) = Vo ~ v( z )i where 77 
is conformal time and z = 1 /a — 1 is the redshift. A ^(k) 
in Eq. (16 1 is the power spectrum of the matter density 


fluctuations today. 

In GR, one usually introduces a fc-independent growth 
function D(z), normalized such that D(z = 0) = 1. The 
density contrast is then given by 6(k, z) = D(z)6(k , 0). In 
addition, due to the absence of anisotropic stress in GR, 
the potentials satisfy $ = , and for subhorizon modes 


In order to calculate Cj 9 and C^ sw using Eqs. ( |16| 
and ( 21 ), we need to compute the time evolution of the 
gravitational potentials 4> and ’E. As discussed in Sec. 
|IV[ we write the equations of motion for the potentials as 
two second-order differential equations (see Refs. [55} [BU. 
for details) and solve them numerically. We set the initial 
conditions at the time of decoupling and then evolve the 
fields to the present time. We additionally assume <I> = 4/ 
at the time of decoupling. The numerical analysis is 
performed using Mathematica. 


VI. COMPARISON TO DATA 

A recent analysis of the ISW cross-correlation has 
been performed in Ref. [96] using the 9-year Wilkin¬ 
son Microwave Anisotropy Probe (WMAP) CMB tem¬ 
perature data [106] . as well as a sample of galaxies and 
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quasars provided by the Wide-field Infrared Survey Ex¬ 
plorer (WISE) 11071 . WISE scanned the entire sky in four 
frequency bands and provided a catalog which contains 
more than 500 million sources, including galaxies with 
redshifts up to z ~ 1 and a median redshift of 0.3 [Ml; 
the dataset is therefore much deeper than the ones pro¬ 
vided by other experiments at similar frequencies. It is 
argued in Ref. [IQ that the large number of sources WISE 
detected, together with the large area of the survey and 
its redshift distribution, makes the catalog one of the best 
for measuring the ISW effect through its cross-correlation 
with the CMB temperature anisotropies. WISE had been 
used previously for the study of the ISW effect (see, e.g., 
Refs. [88] ;HU), but Ref. [96] uses a larger sample by ap¬ 
plying less-conservative cuts to the dataset and a higher 
median redshift. For these reasons, we use the results of 
Ref. [96] in order to compare the ISW predictions of IBB 
to the observations. 

Before moving on to the discussion of our numerical 
analysis and the results, we note that care must be taken 
regarding the bias factor b(z) used in Eq. ( p~7[ ) . Following 
Ref. [96], we use the measured cross-power of galaxies 
and AGN from the WISE survey with the CMB tempera¬ 
ture data. Galaxies and quasars trace the gravitational 
potential up to the bias factor. It is therefore crucial to 
properly measure the bias when comparing the predictions 
of the ISW effect in a theoretical model to the observa¬ 
tional measurements. It is argued in Ref. [56] that the 
best way to measure the bias is through cross-correlating 
lensing potential maps with the tracer field. The reason 
is that gravitational lensing is directly sourced by the 
gravitational potential, and therefore the measurement 
of the bias is more robust and less influenced by system¬ 
atic errors which are present in measurements using the 
auto-correlation spectrum; these errors can lead to an 
incorrect estimation of the amplitude of the ISW effect. 
The authors therefore used weak lensing of the CMB to 
measure the effective bias. The data used for this pro¬ 
cedure are the Planck lensing potential [ 109] as well as 
WISE maps, and the cross-correlation signal is measured 
for 100 < I < 400 (see Refs. [10911110] for earlier studies of 
the correlation between WISE and the Planck lensing po¬ 
tential). Ref. [96] assumes two models for the galaxy bias 
factor as a function of redshift; for a redshift-independent 
bias, b G , their estimated value is b G = 1.41 ±0.15, and for 
a bias linearly dependent on redshift, b G (z) = (1 + z), 

they find b G = 0.98 ± 0.10. For AGN, the bias model is 
b A ( z ) = &q [0.53 + 0.289(1 + z) 2 } with b$ = 1.26 ± 0.23. 
We use the same estimated bias factors in the analysis of 
IBB. 

With the measured values of the bias factor discussed 
above, the ISW effect has been detected in Ref. [96] at the 
level of 3er, with a cross-correlation amplitude, normalized 
to the expected correlation in a ACDM universe, of A = 
1.24 ± 0.47 for the constant bias and A = 1.54 ± 0.59 for 
the linear bias. In both cases 50 million galaxies (with 
a redshift distribution peak at z ~ 0.3) have been used. 
The cross-correlation amplitude using one million AGN 
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FIG. 2. Top panel: Theoretical cross-correlation power 
spectrum Cf 9 between the CMB temperature anisotropies 
and galaxies computed for ACDM and IBB using the density 
contrast today to normalize the weight functions, as well as 
the data points measured using the WISE galaxy catalog in 
Ref. [96]. Bottom panel: The same for the cross-correlation 
spectrum between the CMB temperature anisotropies and 
AGN. 


(with a redshift distribution peak at z ~ 1 . 1 ) has been 
found to be A = 0.88 ± 0.74. 

We now apply the method of calculating the ISW spec¬ 
tra discussed in Sec. [V]to our bigravity model, IBB. This 
is similar to the method employed in Ref. [96] for ACDM. 
We use the same bias factors for galaxies and AGN as 
in Ref. (55], and the redshift distributions dN/dz for the 
galaxies and AGN are taken from Refs. [108 : 1 10 . 

In Fig. [2] we plot the theoretical cross-correlation power 
spectrum Cj 9 between the CMB temperature anisotropies 
and galaxies/AGN computed for both ACDM (fl m = 
0.3) and IBB (/3i = 0.48 and ,84 = 0.94), together with 
the measured values using the WISE galaxy and AGN 
catalogs and the corresponding uncertainties as reported 
in Ref. [55. Here we have used the density contrast today 
to normalize the weight functions ( fl7| an d (|l 8 l ), assumed 
5h = 7.5 x 10 ~ 5 and n = 1 [see Eq. (|20pj, and used 
the constant bias for galaxies. The figure shows that 
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FIG. 3. Theoretical ISW auto-correlation power spectrum 
Cf w for both ACDM and IBB. 


IBB predicts a higher correlation (by roughly a factor 
of 1.5) compared to ACDM, but the amplitudes are still 
consistent with the data. A similar conclusion can be 
made when the weight functions are normalized using 
the gravitational potentials at the time of decoupling or 
when the linearly evolving bias for the galaxies is used. 
Note that our theoretical curves for ACDM are slightly 
different from the ones presented in Ref. |55|, which is 
due to the difference in the choice of transfer function. 
In principle, one could perform the same calculation by 
using parameter values closer to the ones from WMAP9 
or Planck and by choosing a different transfer function, 
but this would only change the theoretical curves slightly. 
Our main objective in this paper is comparing ACDM 
and IBB in terms of the ISW predictions and what is 
important in our analysis is to calculate the theoretical 
curves using the same initial conditions for the two models, 
and infer the relative shift in the ISW amplitudes. 

It is also interesting to see the theoretical ISW auto¬ 
correlation power spectrum C] sw [defined in Eq. ( |2l| )] 
for both ACDM and IBB; this is shown in Fig. [3| There 
are no data points on that plot, since, as we mentioned 
earlier, the signal is too small to be detected directly in 
the CMB spectrum. The figure once again shows that 
the ISW effect in IBB is larger than that in ACDM. It 
is however still subdominant compared to the primary 
CMB signal and is consistent with the data given the 
large uncertainties on large scales due to cosmic variance. 

Before we end this section, let us present another way 
of comparing the ISW amplitudes in ACDM and IBB, this 
time in terms of how well the amplitudes can be measured 
from the same sets of observations. This is based on the 
expected values of signal-to-noise ratio in the two models. 
Following a Fisher matrix analysis presented in Ref. J)(ij . 
the signal-to-noise ratio is: 


/sky ^ Cf T C 9 j + [Cj 9 ] 2 

/sky Ype + 1) T l T , (22) 

where / s k y is the sky coverage fraction, Cj 9 is again the 
cross-correlation power spectrum between the CMB tem¬ 
perature anisotropies and large-scale structure (galaxies or 
AGN), and Cj T and C 99 are the auto-correlation power 
spectra for the CMB and large-scale structure, respec¬ 
tively. The first thing one can infer from this expression 
for the signal-to-noise is that almost all of the contribu¬ 
tions to the ISW effect are from low multipoles and low 
redshifts: they are negligible for l > 100 and z > 1.5 ma¬ 
in Ref. [SB] a theoretical maximum value for the signal-to- 
noise ratio has been evaluated to be 7.6 y/ / s k y for ACDM; 
this is the maximum signal-to-noise one can achieve tak¬ 
ing into account the uncertainties coming from cosmic 
variance. 

Repeating the same calculation of the signal-to- 
noise ratio for IBB, we can compare that value to 
the one for ACDM. We find that the ratio of the 
signal-to-noise as computed for IBB and ACDM, i.e. 
(S'/7V)ibb/(<5'/AI)acdm, is ~1.2 for the galaxies, and ~1.6 
for the AGN. These ratios are closely related to ratio 
of the cross-correlation amplitudes Mibb and Macdm for 
IBB and ACDM, respectively, and show that the differ¬ 
ence between the models is well within the observational 
uncertainties. 


VII. DISCUSSION AND CONCLUSIONS 

In this paper we have studied the integrated Sachs- 
Wolfe (ISW) effect in the Hassan-Rosen theory of bimet¬ 
ric massive gravity. The ISW effect is sensitive to the 
time evolution of the gravitational potentials sourced by 
large-scale structure and is nonzero at late times when 
the cosmic expansion is accelerating. It therefore offers 
a complementary tool for probing the physics governing 
the late-time evolution of the Universe and distinguish¬ 
ing between possible mechanisms generating the cosmic 
acceleration; this includes both dark energy and modified 
gravity theories. The ISW effect contributes to large scales 
(low multipoles) of the CMB temperature auto-correlation 
power spectrum as well as the cross-correlation power 
spectrum between the CMB and tracers of large-scale 
structure such as galaxies and quasars. The contribution 
to the auto-correlation is subdominant compared to the 
primordial contributions and is therefore too small to be 
detected observationally. The cross-correlation contribu¬ 
tions are however large enough and have been detected 
by various studies of the effect. 

We have focused on the infinite-branch bigravity (IBB) 
model, which has viable cosmological solutions in terms 














of both expansion histories and linear growth of structure. 
The model is also free from gradient instabilities and 
can explain the late-time expansion of the Universe with¬ 
out resorting to an explicit cosmological constant, and 
therefore offers an interesting alternative to the standard 
cosmological ACDM model. 

The background and perturbative equations in bigravity 
in general and IBB in particular were reviewed, empha¬ 
sizing the quantities and equations required for our ISW 
analysis. We then presented our method for studying the 
ISW effect in bigravity in terms of the cross-correlation be¬ 
tween the CMB temperature anisotropies and large-scale 
structure, which we employed to calculate the theoret¬ 
ical predictions for a best-fit IBB model and compare 
them to both the predictions in ACDM and observational 
measurements of the effect. 

The data we have used in our analysis are taken from 
Ref. [96], where the 9-year WMAP CMB temperature 
maps as well as the datasets provided by the WISE sur¬ 
vey were analyzed to measure the ISW effect. We used 
the measurements performed in Ref. [96] for the WISE 
catalogs of both galaxies and AGN. We also used the bias 
factors provided by Ref. [55] , 

We applied our method of computing the ISW auto- and 
cross-correlation amplitudes to the IBB model, and com¬ 
pared our results to those of ACDM and to the observa¬ 
tional measurements. We found that IBB predicts higher 
amplitudes compared to ACDM for the ISW temperature 
auto-correlation, as well as the temperature-galaxy and 
temperature-AGN cross-correlations. The amplitudes are 


however consistent with the measurements from WISE, 
and IBB is therefore not ruled out by existing data. 

In order to calculate the cross-correlation power spectra, 
we used two different methods of normalizing the weight 
functions appearing in the expression for the spectra: the 
normalization is done with respect to either the gravita¬ 
tional potentials at the time of decoupling or the matter 
density contrast today. The theoretical amplitudes for 
IBB depend weakly on the choice of the normalization 
scheme, and are in both cases consistent with observa¬ 
tional measurements. We emphasize that ultimately a 
more careful analysis of the evolution of the perturbations 
in a multi-component framework is needed in order to 
set up the correct normalization scheme. This is beyond 
the scope of the present paper and we leave it for future 
work. 
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